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Abstract 

The B® —7- J/ipK® branching fraction is measured in a data sample corresponding 
to 0.41 fb _1 of integrated luminosity collected with the LHCb detector at the LHC. 
This channel is sensitive to the penguin contributions affecting the sin 2/3 measure- 
ment from B° — > J/ipK®. The time-integrated branching fraction is measured to be 
B{B® —7- J/ipK®) = (1.83 ± 0.28) x 10~ 5 . This is the most precise measurement to 
date. 
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1 Introduction 



In the Standard Model (SM) CP violation arises through a single phase in the quark 
mixing matrix pp. In decays of neutral B mesons to a final state which is accessible 
to both B and B, the interference between the amplitude for the direct decay and the 
amplitude for decay via oscillation leads to a time- dependent CP-violating asymmetry 
between the decay time distributions of the two mesons. The mode B° — > J/i/jK® allows 
for the measurement of such an asymmetry, which is parametrised by the B°-B° mixing 
phase (pd- In the SM this phase is equal to 2(3 [2J, where (3 is one of the angles of 
the unitarity triangle of the mixing matrix. This phase is already measured by the B 
factories [3] but an improved measurement is necessary to resolve conclusively the present 
tension in the unitarity triangle fits jl] and determine possible small deviations from the 
SM value. To achieve the required precision, knowledge of the doubly Cabibbo-suppressed 
higher order perturbative corrections, known as penguin diagrams, becomes mandatory. 
The contributions of these penguin diagrams are difficult to calculate reliably and therefore 
need to be extracted directly from experimentally accessible observables. Due to SU(3) 
flavour symmetry, these penguin diagrams can be studied in other decay modes where 
they are not suppressed relative to the tree level diagram. The B° s ->• J/i/>K° mode is the 
most promising candidate from the theoretical perspective since it is related to the I? — > 
J/ipKg mode through the interchange of all d and s quarks ([/-spin symmetry, a subgroup 
of SU(3)) [5] and there is a one-to-one correspondence between all decay topologies in 
these modes, as illustrated in Fig. [T] A further discussion regarding the theory of this 
decay and its potential at LHCb is given in Ref. [6]. 

To extract the parameters related to penguin contributions in these decays, a time- 
dependent CP violation study of the B° s — > J/ip K® mode is required. The measurement 
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Figure 1: Decay topologies contributing to the B° — > J/ipK® and B° s — > J/ipK® channel: 
tree diagram to the left and penguin diagram to the right. 
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of its branching fraction is an important first step, allowing to test the [/-spin symmetry 
assumption that lies at the basis of the proposed approach. The CDF collaboration 
reported the first observation of the — » J/ipK® decay |7j. This letter presents a more 
precise measurement of this branching fraction at the LHCb experiment. 

The strategy of the analysis is to measure the ratio of B° s — > J/ip K® and B° — > J/ip K® 
event yields, which is then converted into a B® — > J/ipK® branching fraction. We make 
use of the B° — > J/ifj K° branching fraction and of the ratio of B° s to B° meson production 
at the LHC, denoted f s /f d [8]. 

We use an integrated luminosity of 0.41 fb _1 of pp collision data recorded at a centre- 
of-mass energy of 7TeV during 2010 and the first half of 2011. The detector [9] is a single- 
arm spectrometer designed to study particles containing b or c quarks. It includes a high 
precision tracking system consisting of a silicon-strip vertex detector surrounding the pp 
interaction region, a large-area silicon-strip detector located upstream of a dipole magnet 
with a bending power of about 4Tm, and three stations of silicon-strip detectors and 
straw drift-tubes placed downstream. The combined tracking system has a momentum 
resolution Ap/p that varies from 0.4% at 5 GeV/c to 0.6% at 100 GeV/c, and an impact pa- 
rameter resolution of 20 [im for tracks with high transverse momentum. Charged hadrons 
are identified using two ring-imaging Cherenkov (RICH) detectors. Muons are identi- 
fied by a muon system composed of alternating layers of iron and multiwire proportional 
chambers. 

The signal simulation sample used for this analysis was generated using the Pythia 6.4 
generator [TU] configured with the parameters detailed in Ref. The EvtGen [T2] . 

Photos [32] and Geant4 [TJ] packages were used to decay unstable particles, generate 
QED radiative corrections and simulate interactions in the detector, respectively. 



2 Data samples and selection 

We search for B — > J/ipK® decays^ where J/0 — » /U + /i _ and K® — > tt + it~ . Events are 
selected by a trigger system consisting of a hardware trigger, which requires muon or 
hadron candidates with high transverse momentum with respect to the beam direction, 
Pt, followed by a two stage software trigger [15j. In the first stage a simplified event 
reconstruction is applied. Events are required to have either two oppositely charged 
muons with combined mass above 2.7 GeV/c 2 , or at least one muon or one high-px track 
(p-T > 1.8 GeV/c) with a large impact parameter with respect to any primary vertex. In 
the second stage a full event reconstruction is performed and only events containing J/ip 
—> candidates are retained. 

In order to reduce the data to a manageable level, very loose requirements are applied 
to suppress background while keeping the signal efficiency high. J/ip candidates are 
created from pairs of oppositely charged muons that have a common vertex and a mass 
in the range 3030-3150 MeV/c 2 . The latter corresponds to about eight times the 
mass resolution at the J/ip mass and covers part of the J/ip radiative tail. The K® 

1 B stands for B° or B°. 
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selection requires two oppositely charged particles reconstructed in the tracking stations 
on either side of the magnet, both with hits in the vertex detector (long K® candidate) 
or not (downstream K® candidate). The K® candidates must be made of tracks forming 
a common vertex and have a mass within eight standard deviations of the K® mass and 
must not be compatible with the A mass under the mass hypothesis that one of the two 
tracks is a proton and the other a pion. 

We select B candidates from combinations of J/0 and K® candidates with mass mj^ K o 
in the range 5200-5500 MeV/c 2 . The latter is computed with the masses of the 
and 7r + 7r~ pairs constrained to the J/0 and K® masses, respectively. The mass and decay 
time of the B are obtained from a decay chain fit [TO] that in addition constrains the 
B candidate to originate from the primary vertex. The \ 2 °f the fit, which has eight 
degrees of freedom, is required to be less than 128 and the estimated uncertainty on the 
B mass must not exceed 30 MeV/c 2 . B candidates are required to have a decay time 
larger than 0.2 ps and K® candidates to have a flight distance larger than five times its 
uncertainty. The offline selected signal candidate is required to be that used for the trigger 
decision at both software trigger stages. About 1% of the selected events have several 
candidates sharing some final state particles. In such cases one candidate per event is 
selected randomly. 

3 Measurement of event yields 

Following the selection described above, a neural network (NN) classifier jTTj is used to 
further discriminate between signal and background. The NN is trained entirely on data, 
using samples that are independent of those used to make the measurements. The training 
maximises the separation of signal and background events using weights determined by 
the sPlot technique [18J. We use the B° — > J/0 K® signal in the data as a proxy for the 
Bg — > J/ip decay. The background events are taken from mass sidebands in the region 
5390-5500 MeV/c 2 , thus avoiding the B° B signal region. A normalisation sample of one 
quarter of the candidates randomly selected is left out in the NN training to allow an 
unbiased measurement of the B° yield. 

We perform an unbinned maximum likelihood fit to the mass distribution of the se- 
lected candidates, shown in Fig. [2j and use it to assign background and signal weights to 
each candidate. The probability density function (PDF) is defined as the sum of a B° 
signal component, a combinatorial background and a small contribution from partially 
reconstructed B — > J/0 K® X decays at masses below the B° mass. The mass lineshape of 
the B° — > J/ip signal in both data and simulation exhibits non-Gaussian tails on both 
sides of the signal peak due to detector resolutions depending on angular distributions in 
the decay. We model the signal shape by an empirical model composed of two Crystal 
Ball (CB) functions [19], one of which has the tail extending to high masses. The two 
CB components are constrained to have the same peak and width, which are allowed to 
vary in the fit. The parameters describing the CB tails are taken from B + — > J/i/jK + 
events which exhibit the same behaviour as B — >■ J/ipK®. The combinatorial background 
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Figure 2: Mass distribution of the B — > J/if)K° candidates used to determine the PDF. 
The solid line is the total PDF composed of the B° — > J/ipK® signal shown in grey and 
the combinatorial background represented by the dotted line. 



is described by a second order polynomial. The B° s — > J/ipK® signal is not included in 
this fit. We extract (14.4 ± 0.2) x 10 3 B° events from the fit. 

The NN uses information about the candidate kinematics, vertex and track quality, 
impact parameter, particle identification information from the RICH and muon detec- 
tors, as well as global event properties like track and primary vertex multiplicities. The 
variables that are used in the NN are chosen not to induce a correlation with the mass 
distribution. This was verified using simulated events. 

To maximise the separation power, a first NN classifier using only the five most dis- 
criminating variables is used to remove 80% of the background events while keeping 95% 
of the B° signal. These variables are the x 2 °f the decay chain fit, the angle between the 
B momentum and the vector from the primary vertex to the decay vertex, the pt of the 
Kg, the estimated uncertainty on the B mass and the impact parameter \ 2 of the J/ip. 

The weighting procedure is then repeated on the remaining candidates and a second 
NN classifier containing 31 variables is trained. A cut is then made on the second NN 
output in order to optimise the expected sensitivity to the B° s yield [20] . 

For the candidates passing the NN requirement, we determine the ratio of B® and B° 
yields for candidates containing a downstream K® or a long K® separately. The B° yield 
is measured in an unbinned likelihood fit to the normalisation sample and scaled to the 
full sample. The B® yield is fitted on the full sample. In both fits, the PDF is identical to 
that used to determine the sWeights with the addition of a PDF for the B® component, 
which is constrained to have the same shape as the B° PDF, shifted by the measured 
B° s — B° mass difference [21]. The results of the fits on the full samples are shown in 
Fig. [3] separately for candidates with downstream and long K®. 
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Figure 3: Fit to full sample after the optimal NN cut has been applied with downstream 
left and long K® to the right. 
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Table 1: B° and B° s yields. Only statistical errors are quoted. The B° yield is obtained in a 
fit to one quarter of the events which have not been used in the NN training (normalisation 
sample) and then scaled to the full sample. 





downstream K® 


long K® 


B° in normalisation sample 


1502 ± 39 


970 ± 31 


B° in normalisation sample (scaled to full) 


6007 ± 157 


3879 ± 124 


B° s in full sample 


72 ± 11 


44 ± 8 


Ratio of B° s to B° 


0.0120 ± 0.0018 


0.0112 ± 0.0020 


Ratio of B° s to B° (weighted average, r) 


0.0117 ±0.0014 



The fitted yields are listed in Table [TJ The long and downstream results are compatible 
with each other and are combined using a weighted average. 

4 Corrections and systematic uncertainties 

Differences in the total selection efficiencies between the B° — > J/ip K® and B® — > J/if> K® 
arise because of the slight difference in momentum spectra of the B mesons and/or the 
final state particles. We find, using simulated events, that the geometrical acceptance of 
the LHCb detector is lower for the B® mode by (1.3 ± 0.5)% where the error is due to the 
limited sample of simulated events. We correct for the ratio of acceptances and assign a 
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Table 2: Summary of corrections and systematic uncertainties on the ratio of branching 
fractions. 



Source 


Correction factor 


Geometrical acceptance (e geom ) 
Trigger and reconstruction 
Decay time acceptance (e t i me ) 
Mass shape 
B° s -B° mass difference 


0.987 ±0.018 
1.000 ±0.010 
0.975 ±0.007 
1.000 ±0.050 
1.000 ±0.004 


Total 


0.962 ±0.053 



conservative systematic uncertainty of 1.8%, which is the sum of the measured difference 
and its error. 

The trigger, reconstruction and selection efficiencies also depend on the transverse 
momentum of the final state particles. Applying the trigger transverse momentum cuts on 
simulated B° and B® decays we find differences of up to 1%, which is taken as systematic 
uncertainty. 

Due to the selection cuts and the correlation of the neural network with the decay 
time, a decay time acceptance function results in different selection efficiencies for the 
B° s and the B°. We determine the lifetime acceptance of the whole selection chain using 
simulated events, and find that the ratio of the time-integrated decay time distributions 
for B° and B° s is 0.975 ± 0.007. The uncertainties on the parametrisation of the lifetime 
acceptance cancel almost perfectly in the ratio, while the ones related to the B° and B® 
lifetimes and the B® decay width difference Ar s do not. 

The largest systematic uncertainty comes from the assumed mass PDF, in particular 
the fraction of the positive tail of the B° extending below the B° s signal. We have studied 
the magnitude of this effect by leaving both tails of the CB shapes free in the fit, or 
by allowing the two CB shapes to have different widths. The maximal deviation we 
observe in the ratios of downstream or long candidates is 5%, which we take as systematic 
uncertainty. The effect of the uncertainty on the B° s -B° mass difference is found to be 
0.4%. 

The corrections and systematic uncertainties affecting the branching fraction ratio are 
listed in Table [2] The total uncertainty is obtained by adding all the uncertainties in 
quadrature. 

We verify that the global event variable distributions, like the number of primary 
vertices and the hit multiplicities, are the same for B° and B° s initial states using the 
B° s — > J/ip cf) channel. We verify that the NN classifier is stable even when variables are 
removed from the training. We search for peaking backgrounds in simulated b — > J/ipX 
events, and in data by inverting the A veto and the flight distance cut. No evidence of 
peaking backgrounds is found. All these tests give results compatible with the measured 
ratio though with a larger statistical uncertainty. 
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5 Determination of branching fraction 

Using the measured ratio r = 0.0117 ± 0.0014 of B° s -»■ J/ipK$ and B° -»■ J/V>#s 
yields, the geometrical (e geom ) and lifetime (e t i me ) acceptance ratios, and assuming 
fs/fd = 0.267 t ao2o [8] we measure the ratio of branching fractions 

£(5° -> J/^ s °) ~ ge ° m tme / s [ ' 

= 0.0420 ± 0.0049 (stat) ± 0.0023 (syst) ± 0.0033 (f s /f d ) 

where the quoted uncertainties are statistical, systematic, and due to the uncertainly in 
f s /f d , respectively. Using the B° -»■ J/^K° branching fraction of (8.71 ±0.32) x 10" 4 [22], 
we determine the time-integrated B° s — > J/ipK® branching fraction 

B(B° S -> J/^Kl) = [1.83 ± 0.21 (stat) ± 0.10 (syst) ± 0.14 (f s /f d ) 

± 0.07 (£(5° J/ipK ))] x 10~ 5 

where the last uncertainty comes from the B° — > J/ipK° branching fraction. This result 
is compatible with, and more precise than, the previous measurement [7j. 



6 Comparison with SU(3) expectations 

It was pointed out in Ref. [23] that because of the sizable decay width difference between 
the heavy and light eigenstates of the B° s system, there is an ambiguity in the definition 
of the branching fractions of B® decays. Due to B° s mixing, a branching fraction defined 
as the ratio of the time integrated number of B° s decays to a final state and the total 
number of B s mesons, is not equal to the CP-average of the decay rates in the flavour 
eigenstate basis 

B(B° S /) theo = T -f {T(B° S -+f) + T(B° S /)) | i=Q , (2) 

used in the theoretical predictions; the restriction to t = removes the effects due to the 
non-zero B s decay width. To obtain the latter quantity from the time-integrated decay 
rates the following correction factor 



A 



J/ipK° 

Ar 



s. 



0.936 ±0.015, (3) 



is applied, where y s = Ar s /2r s is the normalised decay width difference between the light 

and heavy states and A~^ Ks is the final-state dependent asymmetry of the B° s decay rates 
to the J/ipK® final state. In calculating this correction factor we use y s = 0.075±0.010 

and the SM expectation A^}** = 0.84 ± 0.18 
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With this correction, and assuming B(B® — > J/ipK®) theo = \B{B° S — > J/ip K°) theo we 
get the B° s — > J/ipK branching fraction at t — 

B{B° S ->• J/4>K°) thco = (3.42 ± 0.40 (stat) ± 0.19 (syst) ± 0.27 (f s /f d ) 

±0.13(5(5° J/^K )) ±0.05 (y s ,^ArJ) ■ 10~ 5 . 

This branching fraction can be compared to theoretical expectations from SU(3) symme- 
try, which implies an equality of the B® — > J/tpK° and B° — > J/tpn decay widths [6] 

= gCgg -» J/^K°) thco r B o [m B o$(£?° -» J/^tt )] 3 soft) 
where the factor two is associated with the wave function of the 7r°, r R o is the mean B9 s 

D (s) ( s ) 

lifetime and $ refers to the two-body phase-space factors; see e.g. Ref. [5]. 

Taking the measured B(B® — > J/ip i^°) t heo and using the world average [221 EI] for all 
other quantities, this ratio becomes 

3 SE r(3) =0.98 ±0.18 
and is consistent with theoretical expectation of unity under SU(3) symmetry. 



7 Conclusion 

The branching fraction of the Cabibbo-suppressed decay B® — > J/ipK® is measured in 
a 0.41 fb _1 data sample collected with the LHCb detector. We determine the ratio of 
the B° s J/ijKl and B° J/^pKl branching fractions to be fS^#3] = 0.0420 ± 

0.0049 (stat) ± 0.0023 (syst) ± 0.0033 (/,//*)■ Using the world-average B° J/ipK° 
branching fraction we get the time-integrated branching fraction B(B° S — > J/ipK®) = 
[1.83 ±0.21 (stat) ±0.10 (syst) ±0.14 (f s /f d )± 0.07 (B(B° J/ipK ))} x 10~ 5 . The total 
uncertainty of 16% is dominated by the statistical uncertainty. This branching fraction 
is compatible with expectations from SU(3). 

With larger data samples, a time dependent CP-violation measurement of this decay 
will be possible, allowing the experimental determination of the penguin contributions to 
the sin 2/3 measurement from B° — > J/tpK®. 
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